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I. INTRODUCTION 

The presence of quenched disorder in condensed mat- 
ter systems is known to give rise to a wealth of phe- 
nomena, such as new phases and phase transitions, 
new universality classes in critical behavior, localization, 
slow and glassy dynamics, possibly associated with the 
proliferation of metastable states. On the theoretical 
side, quenched disorder requires the introduction of new 
methodological tools. A main technical difficulty comes 
from the fact that a disordered system is intrinsically in- 
homogeneous. As a result, its long-distance properties 
are described by the whole distribution function of some 
disorder-dependent free-energy functional. With the ex- 
ception of rare low-dimensional models for which spe- 
cific methods can be tailored and of phenomenological 
approaches, such the droplet picture, that rely on a real- 
space, hence inhomogencous, description, most theories 
rather deal with disorder-averaged quantities, e.g. cumu- 
lants of a disorder-dependent function or functional. The 
standard tool to achieve this goal is the "replica method" . 

In its most common use, introduced for the study 
of spin glasses by Edwards and Anderson 1 , the replica 
method amounts to computing the average of the log- 
arithm of the disorder-dependent partition function 
of the system under study by introducing n identi- 
cal copies of the original system and, once the av- 
erage over disorder has been performed, by letting 
the number n of "replicas" go to zero with a suit- 
able analytic continuation. The "replicated" hamilto- 
nian, or action in field-theoretical language, is invari- 
ant under permutations of the replicas, i.e. "replica 
symmetric" ; in the standard implementation, external 
fields (or sources), when added, are taken as equal 
for all replicas so that the full hamiltonian/action in- 
cluding the source terms is also replica symmetric-'. 
This replica trick has been mostly used in conjunction 
with mean-ficld-like models (fully interacting, infinite- 
dimensional systems or variables with an infinite num- 
ber of components) : ' ' 1 and approximations (e.g., the 
gaussian variational method ' 6 >'> lS ). More recently, an- 



other variant of the replica method has been considered, 
in the context of the functional renormalization group 
(FRG) for random clastic manifold and random field 
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this case, the fundamental variables of the replicated sys- 
tem are coupled to sources that are different for all repli- 
cas. The source terms then explicitly break the replica 
symmetry. This approach provides a convenient proce- 
dure to generate the cumulants of the appropriate ran- 
dom functional. 

Given the a priori different nature of the two above 
replica schemes, a crucial question is that of their equiv- 
alence in the physical limit where all sources are equal. 
Perturbatively, it is easy to show order by order that the 
two procedures provide the same answer. However, it has 
been realized that perturbation theories fail at describing 
many of the interesting new phenomena induced by the 
presence of quenched disorder, as illustrated by the fail- 
ure of the replica-symmetric solution for the spin-glass 
phase of the Shcrrington-Kirkpatrick (SK) mean-field 
model 2 or the breakdown of the "dimensional reduction" 
predictions in random elastic manifolds and in the 
random field Ising model in low enough dimension 1 . In 
the mean- field-like approaches, in which all replicas are 
coupled to the same source, the way out of the naive, 
perturbative solution comes from a spontaneous break- 
ing of the replica symmetry in the continuation process 
of taking the number of replicas to zero, as in the Parisi 
solution of the SK spin glass model 2 . In the FRG de- 
scription of random manifold and random field systems, 
in which replica symmetry is explicitly broken from the 
beginning by the sources, evading dimensional reduction 
proceeds through the appearance of a nonanalyticity 
in the field dependence of the running renormalizcd 
cumulants of the random free-energy functional as 
the RG flow approaches a zero-temperature fixed 
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This nonanalyticity shows up in the limit where the 
difference between the sources acting on the replicas go 
to zero. Obviously, the perturbative equivalence between 
the two replica approaches, with and without explicit 
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breaking of the permutation symmetry, becomes ques- 
tionable when either one of the approaches encounter a 
singularity. 

The purpose of the present article is then to inves- 
tigate the relation between spontaneous and explicit 
replica symmetry breaking in the theoretical descrip- 
tion of disordered systems. This topic has already been 
studied in a series of articles on random elastic mani- 
folds, starting with the work of Balents, Bouchaud and 
Mezard J9 . In particular, a thorough analysis has been 
carried out by Le Doussal, Wiese et al 1 '' 18 ' 20 in the 
large TV (infinite codimension) limit of the model. In 
this specific case, the authors have been able to unravel 
the full connection between the spontaneously generated 
replica symmetry broken solution previously obtained by 
Mezard and Parisi 3 and the nonanalytic field dependence 
of the cumulants obtained by combining FRG and ex- 
plicit replica symmetry breaking (RSB). The question 
we address in the following is whether a similar connec- 
tion can be established for a more general class of dis- 
ordered systems and for more general, non mean-field, 
conditions. Such a connection would be especially valu- 
able in the case of the random field O(N) model for 
which the existence of a spontaneous RSB phenomenon 
has been suggested 30 ' 31 ' 32,33 whereas recent FRG studies 
have shown that a cusp singularity in the field depen- 
dence of the cumulants of the renormalized random field 
appear below a critical dimension 25 ' 26 ' 27,28 ' 34 . It would 
be interesting too for spin glasses. Indeed, the nature 
of the low-temperature spin-glass phase as well as the 
existence of a transition in the presence of an applied 
magnetic field are still controversial matters. Combining 
the FRG approach and the explicit RSB scheme could 
therefore shed light on the presence or absence of spon- 
taneous RSB in the system. 

When found in theoretical descriptions of disordered 
models, spontaneous RSB in its most common form in- 
volves a "matrix breaking", i.e. involves a 2-point, 2- 
rcplica function, as in the Parisi solution of the SK 
model 2 . As developed in the theory of spin glasses 
and discussed more recently for the random field O(N) 
model 30 ' 31,32 ' 33 , spontaneous RSB is signaled by an in- 
stability of the replica-symmetric solution which is char- 
acterized by the appearance of a zero eigenvalue in 
an appropriate stability operator. The natural frame- 
work to investigate such phenomena involving 2-point, 
2-rcplica functions is the so-called 2-particle irreducible 
(2PI) formalism 35 ' 36 ' 37 : by introducing sources that cou- 
ple not only to the fundamental fields but also to com- 
posite bilinear fields and then performing a double Lcg- 
endre transform, one obtains a generating functional that 
has for argument both fields (magnetizations) and 2- 
point correlation functions. When coupled to a replica 
approach, this formalism allows one to write the above 
mentioned stability operator as the jacobian of the dou- 
ble Legendre transform ! -'> !iX . (On the other hand, the 
FRG approach to disordered systems is commonly ex- 
pressed within the 1-particle irreducible (1PI) formalism 



in which sources are linearly coupled to the fundamental 
fields; the 1PI formalism can be recovered from the 2PI 
one, which makes the latter a suitable starting point for 
studying the relation between spontaneous and explicit 
RSB.) In this work, we show that a general relation can 
be derived between the instability commonly associated 
with a (continuous) spontaneous RSB and the appear- 
ance of a nonanalytic field dependence in the context of 
the explicit RSB approach. This relation suggests a new 
way to tackle putative spontaneous RSB in disordered 
systems, such as spin glasses, by means of a general FRG 
approach. In addition, we illustrate this relation in the 
context of the random field O(N) model, in which we es- 
tablish explicit expressions at leading and next-to-leading 
orders of the 1/TV expansion. 

The article is organized as follows. In section II, we 
introduce the 2PI formalism in the context of the replica 
method for disordered systems. We next compare in sec- 
tion III the replica symmetric and explicit RSB formal 
solutions for the 2-point, 2-replica correlation functions. 
Section IV contains the main result of this study, a re- 
lation between the stability operators that characterize 
spontaneous RSB on the one hand and breaking of ana- 
lyticity in the explicit RSB formalism on the other hand. 
Section V further expands on the consequences of this 
relation and provides some concluding remarks. In ap- 
pendix A we illustrate the relation between stability op- 
erators in the context of the random field O(N) model 
in the large TV limit and in appendix B we consider the 
case of local self-energies. 



II. 2PI FORMALISM WITH REPLICAS 

The 2PI formalism 35 ' 36 ' 37 involves a functional, an "ef- 
fective action", that depends on two kinds of order pa- 
rameters. The first one is the usual order parameter (j) 
- the magnetization - which is relevant to characterize 
the occurence of a paramagnetic-ferromagnetic transi- 
tion. The second type of order parameter consists of 
2-point correlation functions. To implement this formal- 
ism, one has to introduce external sources that not only 
linearly couple to the fundamental variables (as a mag- 
netic field) but also couple to composite fields formed by 
the product of two fundamental variables. 

We consider a system originally defined on a lattice 
and described, at some coarse-grained scale A , by a con- 
tinuous field x( x )- This field represents a spin variable 
for random magnets, a displacement field parametriz- 
ing a manifold in the case of an elastic manifold pinned 
by disorder, etc. We restrict ourselves to a scalar field 
and for simplicity we write down the equations for a 
one-component field; extension to many components is 
straightforward. To this system is associated an effective 
Hamiltonian or action S[x', h] where h is a symbolic no- 
tation to characterize the presence of quenched disorder 
(dilution, random couplings, random magnetic field or 
anisotropy, random potential, etc). The partition func- 
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tion of the system for a given sample with disorder h is 
given by 



Z h [J,K] = / V X exp - S[ X ;h] + / J(x) x (x) 



x(x)K(x,y)x(y) 



x Jy 



(1) 



where the scalar field x is coupled linearly to a source J 
and quadratically to a source K, and where J" denotes 
J d d x. A random, i.e. disorder-dependent, free-energy 
functional Wh [J, K] can be defined as the logarithm of 
Z h [J,K}. 

Within the replica formalism, one introduces n replicas 
of the original system. We consider the general case in 
which the sources are different for all replicas and pairs of 
replicas. After performing the average over the disorder, 
one obtains 

Zrep[{Ja},{ K ab}] = Yl V X 

•* o=l ^ 

n . 
a=l Jx 

+ I I Xa(x)K ab (x,y)xb(y) 



a,b JxJ y 



(2) 

where Latin letters denote replica indices and SVep[{Xa}] 
can be generically written as 

n . 

Srep[{Xa}} = /Z 

a=l Jx 

(3) 

-^E / / S ^Xa,Xb] +••■ 



a ,b Jx J y 



with Stlxa] = S \xa,h], S 2 [xa,Xb] = S\xb]h]S\xa]h] - 
S[\„;h] S[xb',h], etc; an overbar as usual denotes an av- 
erage over the disorder h. Higher-order terms in the num- 
ber of sums over replicas occur when nongaussian distri- 
butions of the disorder are considered. The specific form 
of Si, S2, etc, of course varies with the kind of model 
- random field, random anistropy, random temperature, 
random manifold, spin glass, etc - considered, but at this 
point, we keep the formalism in its full generality. 

From the partition function of the "replicated" sys- 
tem in Eq. (2) one defines the free-energy functional 
W rep [{J a },{K ab }} = \nZ rep [{J a },{K ab }} whose first 
derivatives generate the usual ferromagnetic order pa- 
rameter, 



SW rep [{J e },{K ef }} 



5J a {x) 



(Xa(x)) = (t>aix), 



(4) 



where ( ) denotes a thermal average with the replicated 
action, and the 2-point, 2-rcplica correlation functions, 

SW rep [{J e },{K ef }} 1 

jtt-i — H— = v(Xa{x)Xb{y)} 

= l;{Gab{x,y) + (j)a(x)<f> b {y)) , 

(5) 

where G ab {x,y) = (Xa{x)xb(y)) - (Xa(%)) (Xb{y)) is the 
connected 2-point Green function (or propagator). The 
latter is also obtained through two derivatives of W rep 
with respect to J a (x) and Jb(y)- Higher-order deriva- 
tives with respect to the sources generate higher-order 
connected Green functions. 

The double Legendre transform with respect to the 
two sources J a and K ab defines an effective action 
r 2 p/[{<M,{G ab }]: 

r 2 p/[{(M,{G Qb }] - 

n „ 

-W rep [{J a },{K ab }]+Y^ / Ja{x)<t>a{x) 
a=l Jx 

+ \z2 I f K ab(x,y)(G ab {x,y) + cj) a (x)(j) b (x)), 

Z a,b=l Jx J y 

(6) 

where we have used Eqs. (4) and (5). Differentiation 
with respect to the field <f> a and to the 2-point, 2-rcplica 
correlation function G a b provides 



sr 2P7 [{0 e },{G e/ }] 



$4>a(x) 



Ja(x) + E / K ab (x,y)Mv) (7) 

b J y 



ar 2PJ [{^ e },{G e/ }] 

5G ab (x,y) 



1 



K a b(x,y) 



(8) 



From now on we consider the case where all sources 
{K ab } are set equal to zero. Eq. (8) provides the "equa- 
tion of motion" or Schwingcr-Dyson equation for the 
fields {Gab} which can be formally solved in terms of 
the fields {0 e }: G ab = G* ab [{cb e }}, with 



ST 2PI [{<i> e },{G e f}] 



SG ab (x,y) 







(9) 



Substituting this result into the expression of 
^2Pi[{4>a},{G ab }] leads to the usual, one-particle 
irreducible (1PI), effective action: 



r 2 P/[{M, {G: h [{<M]}] = r 1P/ [{<M] . (10) 



This set of equations provides the starting point of our 



4 



approach. Indeed, Eq. (9), when evaluated for identical 
sources (or in a more restrictive setting for sources and/or 
fields <j) a all equal to zero), provides a self-consistent equa- 
tion for G a b which in limiting (mean-field) cases or with 
approximations (see below) can be used to study the dis- 
ordered model at hand and the occurence of spontaneous 
RSB. On the other hand, when keeping the sources J a 
acting on the various replicas all different, the reduction 
to the 1PI effective action in Eq. (10) suggest the possibil- 
ity of including the FRG studies in the same framework. 

The 2PI effective action is generally written, up to an 
additive constant, as 37 

r 2 p/[{0a},{G afc }] = S rep [{0 a }] + iTrlnG- 1 + 
+ iTrGG 1 ({<M)+ r 2 [{<M,{G Qb }], 

(11) 

where the trace involves a sum over replica indices and an 
integration over space; Gq 1 is the classical (bare) inverse 
propagator: 

and r2[{</> Q }, {G a {,}] includes all 2PI components of the 
effective action. To see this, one derives Eq. (11) with 
respect to the G ao 's (again at vanishing sources K ao ). 
One finds 

[G*- 1 ({^})] ah = [Go^eDL 

(13) 

-E a6 [{<M;{G c * a ({0 e })}] 

where we have defined 

v r rj. \ in w - r /r 2 [-Q e }, {G cd }} 

Xab[x,y;{j> e },{Gcd}\ = -2 jtt-, r , (14 

oGabix, y) 

which is nothing but the self-energy for the replicated 
system. As well known, the self-energy E is made up of 
1PI diagrams, and a functional derivative with respect 
to G corresponds to opening one propagator line. This 
implies that r 2 is made up of 2PI vacuum diagrams only, 
with dressed propagator G. 

The 2PI functional is of course not known exactly ex- 
cept for very specific conditions. Approximations are 
thus required for practical computations on disordered 
systems. In many approximations used in the literature, 
the self-energy is taken as purely local and all momentum 
dependences are contained in the bare inverse propaga- 
tor Gq . This is verified in mean- field models 2 ''' and is 
the case for (i) the Gaussian variational method which 
amounts to replacing the original Hamiltonian by a vari- 
ational Gaussian one°' fa ' 7 ' 8 (and is exact in the limit of 
large number N of spin components' 5 ) and (ii) for the 
"locator" approximation 39 ' 40 ' 41 ' 42 and related dynamical 
mean- field theory (DMFT) approaches ! : , in which the 
local part of the self-energy is calculated through an ef- 
fective single site model. Other types of approximations 



involve expansions either in loops or in powers of 1 /N in 
the large TV limit. We provide examples in Appendices 
A and B. 



III. REPLICA SYMMETRIC AND EXPLICIT 
RSB SOLUTIONS 

Consider first the situation in which all sources J a are 
equal, J a = J (recall that all sources K a b are taken to 
zero) . The replicated action including the source terms is 
now invariant under permutation of the replicas. When 
no spontaneous breaking of replica symmetry occurs, the 
solution of the (Schwinger-Dyson) equation of motion, 
Eq. (13), is replica symmetric (RS), namely <f> a = </>, Va, 
G a b{<t>] = Gc{4>]5 a b + G £>[<!>], where Gq and Go as usual 
denote the (disorder) connected and disconnected pair 
correlation functions (from now on we drop the star that 
characterizes the solution) . The algebra of RS matrices is 
rather simple (see the above form of G a b) and by further 
restricting our study to a uniform source, Eq. (13) can 
be rewritten as 

G- C \q- 0) - G^(q; 0) - E c [q; 0; G c {4>), Gd (</>)] , (15) 

G D \q; <f>) = G^(q; 0) - E D [q; 0; G c (<t>), G D (cf>)] , (16) 

where we have also decomposed the RS self-energy ac- 
cording to 

E a6 [q; <f>; G c (0) , G D (0)] = E c [q; <t>; G c (</>), G D (</>)} S ab 

+ E D [q-fr G c (<i>),G D (4>)] 
(17) 

and have kept square brackets to indicate that the self- 
energies are still functionals of the correlation func- 
tions Gc, Gjj. The inverse connected and disconnected 
functions (i.e., the connected and disconnected 2-point 
proper vertices) are given, in the limit where the number 
of replicas goes to zero, by 

Gc 1 (q;ct>) = Gc(q;<t>)- 1 , (18) 

G D \q- 4>) = -G c (q; 0y 2 G D (q; 0). (19) 

The solution of the Schwinger-Dyson equation in the 
presence of an explicit RSB due to the sources J a 
is a priori much harder to handle, even at a for- 
mal level. However, a workable scheme is provided 
by using systematic expansions in increasing number 
of unrestricted sums over replicas (or "free replica 
sums") 11 ' 16 ' 17 ' 18 ' 25 ' 26 ' 27 - 28 ' 44 . For instance, any function 
of the replica fields <f) ai say a function F({<f> a }), can be 
expanded as 

F({W)=Ei £ F p(<t>ai,-,K) (2°) 
p>l a 1; ...,a p 
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where the _F p 's are continuous and symmetric func- 
tions of their arguments, whereas matrices, say a matrix 
A ab ({4>e}) can be written as 

A ab ({<P e }) = A a ({<j> e }) S ab + A ab ({<p e }), (21) 

where A ab does not contain any Kroneckcr symbol and 

A a {{<p e }) =J?°]faO 

p>l ' ei ,...,e p 
2 ab ({0 e }) =Al°](0 aj t ) 

p>l ei ,...,e p 

(23) 

Higher-order tensors can be expanded as well . Terms 
appearing in the various equations can then be expanded 
in number of free replica sums and, after properly col- 
lecting and symmetrizing, one can proceed to an order- 
by-order identification. 

Applied to Eq. (13) and again to uniform sources J a , 
the procedure leads at zeroth-order to the following ex- 
pressions: 

G^iq-Aa) = G [ ° ] -\q^ a ) £ [0 W a ), (24) 



the same as those that are arguments of the correlation 
functions entering in the functionals of the right-hand 
side of Eqs. (26,27). This is a direct consequence of the 
expansion in free replica sums. In addition, the functions 
depending on 2-replica fields, G^ (q; cf> a , 4> b ) and the as- 
sociate proper vertex and self-energy, are invariant under 
the exchange of the arguments and are therefore even in 

4>a - <Pb- 

Suppose now that all replica fields {4> a } arc equal, 
4> a = (f> } Vo. Provided that the dependence of the func- 
tions on the replica field arguments is analytic when 
the arguments become equal, any matrix A ab ({(j) e }) is 
RS. As a consequence, it can as well be written as 
Ac(4>)S a b + Ad{4>) and, in the limit where the number of 
replicas goes to zero, one finds that 

A a ({4> e = <{>})= A® (<t>)=A c (4>), (30) 



A ab {{4> e = </>}) = AW{4>,<t>) = A D {4>). (31) 

In particular, G^(q;<f)) = G c (q;4>), &°}{q\4>,<f>) = 
G D (q;4>), and Eqs. (24-29) coincide with Eqs. (15-19). 
This property, as we shall delve on more below, relies on 
the assumption of analytic field dependence in the limit 
where all fields arc equal. If not the case, Eqs. (24-27) are 
a priori more general than the RS equations, Eqs. (15- 
17). 



where 



;0 Q ;GM(0 Q ),G[ o ](0 Q ,0 a ) 



(25) 



(26) 



£[°]( g ;0 Q ,0 o ) = £ 



[o] 



g ;0 a ,0 o ;G[°](0 a ),G[ o ](0 b ), 



G[°](0 a ,^),G[ o ](^,0 6 ),G[ o ](^,</) 6 ) , 

(27) 

and the inverse functions G^ -1 , G^ ' -1 are given by 

G^- 1 (q;cf >a ) = GW(q;cf, a r\ (28) 



G^- l {q^ a Ab) 



G^tq-^G^tq^aAb^WAb) 



_ x (29) 



Note that E^ and E^ are still functionals of the 2-point 
functions G^l and G^ . However, an important observa- 
tion is that the replica fields a , <fi b that are arguments of 
the self-energies in the left-hand side of Eqs. (26,27) are 



IV. STABILITY CONSIDERATIONS 

We now establish the main result of this article, i.e. 
the relation between the existence of an instability of the 
RS solution for the 2-point, 2-replica functions and the 
appearance of a nonanalyticity in the field dependence 
of these functions in the explicit RSB approach. The 
2PI formalism allows one to investigate the stability of 
the solutions to the Schwingcr-Dyson equations, either 
in the RS scheme [Eqs. (15-17)] or in the explicit RSB 
scheme [Eqs. (24-27)]. 



A. Stability of the RS solution 

We start by looking at the stability of the RS solution 
when n — ► 0. This involves considering the eigenvalues 
of the matrix formed by all the second functional deriva- 
tives of T2P1 with respect to the 4> a 's and G ab 's. As is 
well documented since the early work of de Almeida and 
Thouless on the Sherrington-Kirkpatrick mean-field spin 
glass model 45 , an instability with respect to a continuous 
breaking of replica symmetry is signaled by the fact that 
the stability matrix in the so-called "replicon sector" is 
no longer strictly positive definite. Restricting ourselves 
to the study of this replicon component, we only need 
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the matrix operator 



A ab , c d{q,q') 



q" 



S(q' 



(32) 



<5 2 r 2P/ [{«/» e },{G e/ }] 



SG ab (q, q")5G cd {q',q" 



RS 



where the subscript RS indicates that the derivative is 
evaluated at the RS solution of Eqs. (15-17), / denotes 

J d d q/(2n) d , and the replica indices are such that a^b 
and c 7^ d; the factor of 2 is introduced for further 
convenience (in the derivation in the right-hand side of 
Eq. (32), G ab and G ba are first formally considered as dif- 
ferent). With uniform souces J a and vanishing K ab s, the 
system is translationally invariant which brings a global 
factor 5(q + q' + q" + q'" ) in the second derivative of T 2 pi ; 
the delta factor in the right-hand side could therefore be 
replaced by S(q + q") as well. It is convenient to con- 
sider a symmetrized form of A abtCd which, due to replica 
symmetry, has only 3 different components, 



A\(q, q') = sym{A a b,ab(q,q')} 
A%{q,q') = sym{A afc . QC (g,q')} 
&3(q, q') = sym{A ab . cd (q, q')} 



(33) 



where sym{A ab . ab } indicates (l/4)(A ab . ab + A ba , ab + 
A a 6,6a+A6 a ,(,o), etc; A x , A 2 and A 3 respectively involve 2, 
3 and 4 distinct replica indices. The rcplicon eigenvalue 
of the matrix is then given by 

A R (q, q 1 ) = Ai (q, q') - 2A 2 {q, q') + A 3 (q, q'). (34) 

Note that in the rather generic case considered here, the 
replicon is an operator that is momentum dependent. It 
is however not the most general form as we have some- 
what restricted the momentum dependence in Eq. (32). 
This is nonetheless sufficient when the system of interest 
is statistically uniform (which is the case of the replicated 
system described by Eq. (2)) with uniform and transla- 
tionally invariant applied sources . 

The replicon operator is definite positive in the param- 
agnetic phase that is generically present at high enough 
temperature and disorder strength. The vanishing of its 
smallest eigenvalue, if it indeed occurs, defines a gen- 
eralization of the de Almeida-Thoulcss line of instabil- 
ity found in the SK spin-glass model . In the lat- 
ter case, it characterizes the instability of the RS so- 
lution toward spontaneous RSB. Such an interpreta- 
tion also applies to the random elastic manifold model 
in the large N limit (i.e., infinite codimension for the 
manifold) 1 ' ' 18, 44 and more generally to the systems stud- 
ied via either the gaussian variational method 5 ' 8 or the lo- 



cator approximation 



39,40,41,42 



. In the above cases though, 



the mean-field-like character of the problem or of the ap- 
proximation simplifies the momentum dependence: it is 



sufficient to consider local self-energies and (spin glass) 
order parameters so that the replicon operator is easily 
diagonalized (see Appendix B). The replicon eigenvalue 
is then related to the so-called spin-glass susceptibility, 
which is defined in the original disordered model as 



XSG 



1 

V 



((x(x)x(y))-(x(x))(x(y))) 2 (35) 



and is obtained, in principle at least, in measurements on 
a single macroscopic sample. In this case, the spin-glass 
susceptibility diverges when the replicon eigenvalue goes 
to zero (see Appendix B). 

The more general setting considered here also en- 
compasses the random field 0(N) model as previously 
discussed ' . The instability associated with the vanish- 
ing of the smallest eigenvalue of the replicon operator is 
harder to interpret than in mean-ficld-likc cases and does 
not necessarily implies the divergence of the spin-glass 
susceptibility (see below). It has for instance been inter- 
preted as signaling the appearance of bound states be- 
tween pairs of replicas, the replicon operator being then 
akin to the corresponding Bethe-Salpeter kernel 4 ''. 



B. Validity of the analytic solution in the explicit 
RSB formalism 



We now investigate the stability of the solution of 
Eq. (25) with respect to the possible appearance of a 
nonanalytic field dependence when the two field argu- 
ments 4>a and 4> b become equal. To this end we apply the 
operator d^ a — d^ b (with d r p a = d/dcj) a ) to both sides of 
Eq. (25) and use Eqs. (27) and (29). We obtain, under 
the natural assumption that the bare inverse propagator 
Gq ' 1 is regular when <f) a — > <p b : 



G^iq-^G^iq-M-Hiq 



<5£[°](g;^,0 b ) 



(d 0a -^JG[°V;^,</> b ) 



<5G[°V;0 a ,0 b 

(36) 



in which the right-hand side involves derivatives 
of the self-energy S^l with respect to its explicit 
dependence on <p a and <p b as well as with respect 
to d^4 a ),Gl°\<k),GM(4> a ,<f> a ) and &°\^ b ^ b ), 
keeping G'- ' ((f) a , (f> b ) constant, whereas the deriva- 
tive of S[°l in the left-hand side is taken with 
4> a , <f> b , G[°l(0a), GM (0 b ), GM (0 a , (j> a ), Gl°l {<t> h , 4> b ) con- 
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stant. We define the operator A^l as 
kW{q,q';ct> a ,ct> b ) = 
G^{q;<p a )- x G^{q;4> b )-H{q + q')- 



6EW(g;<l> a ,<j> b ) 
JG[°](g';0 o ,0 6 ) 



(37) 



Consider now the limit of Eq. (36) when <j) a = <p b = <fi. 
Due to the symmetry properties (recall that the 
functions are even in tfi a — <f> b ), the right-hand side 
of the equation vanishes, as one does not expect in 
this expression any singular behavior in the field 
arguments. (In doing so, we are guided by the pre- 
vious work on nonanalytic field dependence in disordered 

Systems 10,11 ' 12 > 13 > 14 > 15 > 16,17,18, 19,20,22,23, 24,25,26,27,28,34,44,47 

that show that a singularity is indeed only expected in 
the 4> a — (j>b dependence of &°\q'; <f> a , </>&)). Eq. (36) then 
becomes 

A[°W;0» (fy. -fyJG[°V;<k,&)|* = ' 

(38) 

where the derivative is evaluated for <f> a = <f> b = 4>. 
If A[°l(g, q ; ; <f>, cf>) is a strictly positive definite opera- 
tor, the only solution of the above equation is (<9^ a — 

d < p b )&°'(q f ; <f> a , <t>b)\^ = 0. Similarly, all higher-order odd 

derivatives of G ' with respect to 4> a — 4> b (keeping e.g. 
4>a + 0b constant) vanish, which characterizes an analytic 
dependence of the function on 4> a — <p b near the equality 
of the 2-replica fields. This is what is expected in the 
generic paramagnetic phase found at high enough tem- 
perature and disorder strength. 

On the other hand, if a cusp is present in the field 
dependence of G ' , namely, if 



GE P (g;0)|0 a -0 b | Q + 



(39) 



when (f> a ,(f> b — > 4> with < a < 2, the operator 
A[°](g, q'; (f>, <f>) must have a vanishing eigenvalue so that 
Eq. (38) has a solution with {d^ a —d^ b )&^ (q 1 ; <j) a , <fo)L ^ 
0. The same actually applies to any nonanalytic de- 
pendence on <f> a — <f> b since it implies that a given odd 
derivative of the function (under the action of d$ a — ) 
is nonzero despite the fact that the function is even in 
4>a — <f>b- Conversely, if the operator A[°l(g, q'; (f>, (f>) be- 
comes marginal with a vanishing eigenvalue, the function 
G^ may be written in the vicinity of (j) a = <f) b = <j> as 

Gl°] (g; <fa, <M = CM (g; <f>, 4>) + F(<j> a - <M *o(«; 4>) + 

$±(q;<f>a,<f>b), 

(40) 

where ^o(q', 4 1 ) is the eigenvector associated with the van- 
ishing eigenvalue of A^ (g, q'; <f>, 0), F(<j) a — <j) b ) is an even 



and possibly nonanalytic function of its argument, and 
is orthogonal to ^o{Q',^>), regular in its field argu- 
ments and even in <p a — 4>b ■ It is easy to verify that the 
above expression Eq. (40) satisfies Eq. (38) and all higher- 
order equations involving odd derivatives. The solution 
must of course satisfy Eq. (36), which puts constraints on 
the functional form of F((f> a — <p b ) . One can however won- 
der if the Schwingcr-Dyson equation Eq. (25) combined 
with Eq. (27) is still valid when a cusp in |</> a — 0b| has ap- 
peared in G'°l . This question has been examined in detail 
in the large N version of the random manifold model at 
leading order 1 '- 18 ' 20 where it has been shown that a cusp 
indeed occurs. In this case, the analog of Eqs. (25,27) 
have no solution in the "cuspy" region. However, the 
main point of the analysis of [17,18,20] is that when put 
in a differential form, via a FRG approach, the equations 
can then be continued and solved. 



C. Proof of equivalence 

To prove the identity of the operator A^l (g, g'; 0, <j>) 
signaling the limit of validity of the analytic solution 
in an explicit RSB scheme and the replicon operator 
Afl(g, g'; 4>) of the RS scheme, it is convenient to rederive 
AM by applying the operator d^ a — d^ b to the general 
Schwingcr-Dyson equation, Eqs. (9,13), and by using an 
expansion in free replica sums only at a later stage of the 
process. 

Applying the operator d<p a — d$ b to Eq. (9) with a =^ b 
and with uniform sources leads to 



= {d<j> a -dfo) 



<5r 2P7 [{</» e },{G e/ }] 



8Gab{q) 



E 



<5 2 r 2P/ [{0 e },{G e/ }] 



8G ab {q)8G cd (q') 



(41) 

where we have momentarily reintroduced a star to indi- 
cate that the quantities are evaluated with the solution of 
the Schwinger-Dyson equations (symbolically, G = G*). 

We first study the first term of the right-hand side of 
Eq. (41). From Eqs. (11) and (14) one has 

ST 2PI [{^ e },{G ef }} 
SG ab (q) 

-G- b 1 (q) + G^ b (q;{4> e }) + ^ab(q;{cf>e},{Gef}) ■ 

(42) 

We apply d$ a — d<p b to this expression and evaluate the 
result both at the solution G — G* and in a configura- 
tion where all the replica fields </> e are equal to <p. Under 
this operation the first term in the right-hand member 
of Eq. (42) trivially vanishes since the G e /'s and the e 's 
are considered here as independent variables. The second 
and third terms also vanish under the assumption that 
the dependence of G^ b (q; {4> e }) and S ab (g; {0 e }, {G e /}) 
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on the </> e 's (at constant G e /'s) is regular: the two func- 
tions being symmetric in the exchange of the two fields 
4> a and 4>b, the difference of derivatives vanish in the limit 
where all replica fields become equal. 

We now consider the second term of the right-hand side 
of Eq. (41), which we also evaluate in a configuration with 
all replica fields equal to 4>. Introducing 



M 1 (q,q') = 2 



<5 2 r 2P/ [{0 e },{G e/ }] 



5G ab (q)6G aa (q') 

5 2 T2Pl[{<j>e},{Gef}] 



and, when c ^ a, 6, 



M 2 (q,q')=2 



5G ab (q)5G bb (q>) 



<5 2 r 2P/ [{<M,{G e/ }] 



5G ab (q)SG cc (q') 
we can thus rewrite Eq. (41) as 



(43) 



(44) 



bounded and that no singularities are encountered in the 
free sums over replicas, e.g. that 



a*. E G[0] (9; fa > &)] L = d ^ E G[0] ^ &)] I 



one readily finds that 



G: a (<7;{^}) + G^( g; {0 e }) 



= 0, 



(47) 



(48) 



(<9«/> a - a 06 )G; c (g;{<?!) e })| = 0, whenc^ a, 6, 



so that the first term in the left-hand side of Eq. (45) 
vanishes and the only surviving term in Eq. (46) is the 
one involving the combination Ai — 2A 2 + A3. As a final 
result, Eq. (45) then reduces to Eq. (38) with 

AP% q'; 0, 0) = Ai(q, q'; 0) - 2A 2 (q, q'; 0) + A 3 (<?, q'; 0) 



Mi(q,q') (fy. - d*,) [G* aa {q';m) + G & * b (</; {0 e })] 



+M 2 (q,q') E (^-^)G: C ( 9 ';{«)| 



+ J i (d^-d <f>b ){2A 1 (q,q')G* ab (q';{4, e }) 



-2A 2 ( g ,g') E 



G !c(<l'; WJ) + G* bc (q'; {fa}) 



+A 3 (q,q') E G *c*(?';{<«) 



(45) 

where the A,'s, i = 1,2,3, are evaluated in a configura- 
tion with all replica fields equal (i.e., replica symmetric) 
and are given by Eqs. (32,33). They can then be moved 
under the action of the operator d$ a — d$ b . The sec- 
ond term of the left-hand side of Eq. (45) can be further 
recxpressed as 



2 (A x - 2A 2 + A 3 ) (d^ - d^) G* ab (q'; {cj> e })\, + 



2(A 2 -A 3 ) (d^-d^) 



Eg: c (</;{<M) 



E G* bc {q'; {4> e }) + A 3 (d^ a - d 4b ) E G*M\ We}) 

(46) 

We finally expand the different correlation functions 
G* ab {q;{4>e}) in number of free replica sums. Under 
weak assumptions, namely that all first derivatives are 



(49) 



The A, 's are equivalently obtained as zeroth-order contri- 
butions in the free replica sum expansions for a configu- 
ration of equal replica fields or from the RS solution when 
n — ► (the two approaches being equivalent as discussed 
in section III). Eq. (49) is our main result: it proves 
the equivalence of the operator signaling a breakdown of 
the analytic solution in the explicit RSB scheme and the 
replicon operator associated with the stability of the RS 
solution in the standard replica scheme. The demonstra- 
tion can be easily extended to systems described by an 
A^-component field. 

In Appendix A, we consider the random field O(N) 
model in the 2PI formalism at leading and next-to- 
lcading orders of the 1/N expansion. We explicitly 
show that the operator associated with the instability 
(or breakdown) of the analytic solution in the explicit 
RSB scheme reduces to the replicon operator derived in 
Refs. [31,32] in the appropriate limit. 



V. DISCUSSION 



A. Susceptibilities and spin-glass behavior 

As pointed out in Ref. [32], the inverse of the repli- 
con operator, XB,(q,<l') = A^ {q,q'), plays the role of a 
susceptibility (operator), and a continuous "transition" 
associated with spontaneous RSB is characterized by the 
divergence of the largest eigenvalue of the susceptibility 
operator. Similarly, according to the arguments of sec- 
tion IV, such a divergence also signals the appearance 
of a nonanalytic field dependence of 2-rcplica correlation 
functions. It is instructive to relate this susceptibility op- 
erator to other quantities. On the one hand, XRilyl') i s 
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connected to the standard spin-glass susceptibility, which 
is defined in Eq. (35), as 32 

Xsg= f [ Xn(q,q')- (50) 

J q J q' 



On the other hand, xn(q, ?') can be expressed in terms 
of usual many-point correlation (Green) functions that 
are generated by derivatives of the replicated free-energy 
functional W rep with respect to the sources J Q (x), in the 
absence of any source K a b(x,y) coupled to composite 
replica fields. After introducing 

q', q"; 0) = %' + q") G^ cd (q, q'; 0), (51) 

where WifL is the 4th derivative of the generating func- 
tional W with respect to J ai Jb, J c , Jd and is evaluated 
in a RS configuration with all replica fields equal to 0, 

and 28 

(52) 

G { a l b (q, q'\<t>)~ 2GH(<z, q'; 0) + G a % d {q, q';<f>), 

where distinct Latin indices denote here distinct replicas, 
one finds that 

XR(q, q') = Gf (q, <?';<£) + S(q + q') G(q; <f>) 2 , (53) 

with, we recall, G(q;<fi) = Gc{q]4>) (see section III). If 
one expresses the 2-replica propagator in the explicit RSB 
scheme at zeroth order, G^(g), not as a function of the 
replica fields <p a , <pb but as a function of the (uniform) 
sources J , Jf,, one also derives that 

dj a dj b (q;J a ,J b )\ J = G<£ (q, 0; 0) = Gf (0, q; <j>) 
= XR(q,0)-S(q)G(q;<f>) 2 , 

(54) 

where the partial derivative is evaluated for J a = Jf, = J. 
From this, it can be seen that if the passage from <f> to 
J is nonsingular, i.e. if the ferromagnetic susceptibility 
Xfm = G(q — 0) stays finite, the existence of a cusp in 
4>a — 4>bi or alternatively in J a — Jf,, implies a divergence 
of Xn(q,0) = XR(0,q)- 

This exercise may prove interesting in the case of dis- 
ordered systems for which it can be shown that the spin- 
glass susceptibility is bound from above by the ferromag- 
netic susceptibility, xsg < Xfm, and therefore cannot 
diverge if the latter stays finite. This property has re- 
cently been shown for the random field Ising model . In 
such cases, no spin-glass phase, defined as a phase with 
a diverging xsg, can be found outside a paramagnetic- 
ferromagnetic critical point. However, this argument 



does not preclude the existence of a phase with spon- 
taneous RSB and/or a phase with cuspy 2-replica cor- 
relation functions. What is required for such a phase is 
the divergence of the largest value of the replicon sus- 
ceptibility operator. If this eigenvalue is of measure zero 
in the spectrum, its divergence may not necessarily trig- 
ger that of XSG (see Eq. (50)). On the other hand, it 
could still imply the divergence of x/?(g, 0). This is ad- 
mittedly a quite intricate scenario. Calling "glassy" such 
a phase with anomalous behavior but no diverging XSG 
has anyway little physical motivation unless the glassy 
characteristics are supported by other phenomena such 
as slow dynamics. 



B. Renormalization group 

Formally relating the stability conditions of the RS so- 
lution in the standard replica method and of the analytic 
solution in the explicit RSB scheme is of course only one 
piece of the puzzle. It remains to be seen for specific dis- 
ordered systems if RSB is actually spontaneously broken 
and if a cusp actually appears in the field dependence 
of the 2-replica pair correlation functions. One potential 
advantage of the explicit RSB scheme is that it can be 
combined with renormalization group ideas, more specif- 
ically with the FRG. In disordered systems such as ran- 
dom manifold 9 ' 10 ' 11 ' 12 ' 13 ' 14 ' 15 ' 16 ' 17 ' 18 ' 20 - 21 ' 22 , random field 
and random anisotropy 9 ' 19 ' 23 ' 24 ' 25 ' 26 ' 27 ' 28 ' 34 ' 47 models, it 
has indeed been possible to generate an infinite hier- 
archy of flow equations for the renormalized cumulants 
of the disorder (which are related to the many-replica 
correlation functions for nonequal sources 2 '). On the 
other hand, implementing RG schemes to take into ac- 
count fluctuations and go beyond mean-field treatment 
has proven difficult in the presence of spontaneous RSB 
(see e.g. Ref.[38,49]), which, for instance, still prevents 
one from reaching definite conclusions concerning finite- 
dimensional spin glasses. Looking for cuspy behavior in 
an FRG formalism with explicit RSB therefore appears 
as a promising alternative. 

The presentation given in the previous sections of 
the explicit RSB method within the 2PI formalism has 
however not incorporated any RG scheme. RG can 
be implemented through the introduction of appropri- 
ate cut-off functions that suppress the contribution of 
the low-momentum modes below some running scale 
£50,51,52,53,54,55 _ This, a s well as the connection with the 
1PI FRG formalism developed in Refs. [25,26,27,28,34], 
will be discussed in a forthcoming publication ' 1, . For 
now, we note that quite generally a cusp is associated 
with "shocks" or "avalanches" 21 ' 22 ' 28,29 which correspond 
to abrupt changes in the ground state of the system. 
(We focus here on equilibrium conditions, hence on the 
ground state, not on metastable states and hysteresis 
phenomena"). At zero temperature, for many disor- 
dered systems (random field and anisotropy, spin glass, 
etc.), such abrupt changes are generic (see e.g. the sim- 
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ulation results in Refs. [58,59,60]) so that one indeed ex- 
pects to observe cusp behavior in the field dependence of 
the 2-replica correlation functions. However, important 
questions are then whether such a phenomenon survives 
at long distances, e.g. at the critical point in random- 
field systems, and if it persists at nonzero temperature 
where finite-scale discontinuities are usually believed to 
be rounded. FRG is a tool to answer such questions. <>4 



APPENDIX A: 1/N EXPANSION FOR THE 
RANDOM FIELD O(N) MODEL 

The 1 /N expansion of the 2PI effective action has been 
developed mainly in the context of ouf-of-cquilibrium dy- 
namics to derive nonperturbative approximations. The 
rules defining this expansion have been dicussed at length 
in Ref. [37,61,62]. We only recall here the main features 
of this expansion and we closely follows the notations 
employed in Ref. [61,62]. 

The action is assumed to be of order N and the fields 
of order y/N. The 2PI effective action is a singlet under 
the rotation group O(N). It must be constructed from 
the O(N) invariants: (/> 2 , TrG p and Tr^ 2 G p . We write 
the 1/N expansion of I^^G], which is the sum of 2PI 
diagrams in the expression of Y2PI (see Eq. (11)), as 



T 2 [<t>,G\ = T^[<t>,G\+T» LU \<t>,G\ 



(Al) 



where T^°[<f>,G\ and T^ LO [<p,G] correspond to the lead- 
ing order and next to leading contributions that gather 
terms of order N and 1, respectively. 

We start with the replicated action for the random field 
O(N) model 30 : 



S[{Xa}} = [ WilUdXa) 

Jx ^ a=l V - 



m 2 \Xa\ 2 



A 



12AT 



\Xa\ 



The leading-order (LO) contribution to T2[(/>,G] is sim- 
ply obtained from its counterpart in the O(N) case, up 
to the introduction of replica indices: 

r*°[*,G] = 4^E/ (S^fr*)) • ( A4 ) 



The above expression has a graphical interpretation 
which is a direct extension of what is used in the stan- 
dard O(N) model. Although graphs provide a compact 
representation, we shall not give them in this appendix. 

The next-to-leading order (NLO) is also given by a gen- 
eralization of the O(N) case. To derive it, it is convenient 
to consider separately the ^-independent part r^nd an d 
the (^-dependent one r^ 1 ^ . The former is generated by 
the vertex associated with XaXaXaXa- This contribution 
can be formally written as 



IS[G] = ^Trln(B) 



(A5) 



with 



B ab {x, y;G) = S(x - y)8 ab + — ^ &%(x, y) G%(x, y) 

(A6) 

and where the trace acts both on replica indices and spa- 
tial coordinates. 

In turn, the ^-dependent part of the NLO contribu- 
tion is generated by the vertex associated with (f>aXaX b X b 
with 4>£ being considered as an inserted operator. It is 
given by 



r^ e °>,G] = 



y - 

^ 2 

a,b=l 



Xa -Xb 



(A2) 

where \a is an A^-componcnt replica field with a = l...n 
and A is the bare second cumulant of the (Gaussian dis- 
tributed) random field. Latin indices are used for replicas 
and Greek indices for vector components. 

For uniform replica field configurations a , the free 
inverse propagator is then given in momentum space by 



6^ a|2 



+—J>XSab-A5^}5( q + q ') 



= [Go Tabil) <% + <?'), 

where [Go]^ are the components of Go, a &- 



(A3) 



/ / ^U(x,y;G) y £<%(x)G%(x i v)<P b (v) 



x J y a,b 



I'- 1 ' 



(A7) 



with l ab (x,y;G) iteratively defined as 



l ab (x, y; G) = — ]T G%{x, y) G%(x, y) 



ft J J 



/E lac(x >Z ;G)J2G^(z,y)G^(z,y) . 

(A8) 

Note that l a b(x,y;G) and B a b(x,y;G) are related by 
B-t(x,y;G) = 6(x - y)6 ab -I ab (x,y;G). (A9) 
Finally, gathering all LO and NLO contributions and 
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considering uniform replica field configurations leads to 
V2&.G A 



V 4\N 



E(/E g m) 

4E/ Wfl)L(?) 

EE« / W?;G) G£(g). 

a,b 9 



(A10) 



It is convenient to introduce a "polarization" 



IU(,;G) = A / E G ^(9')G^Cz + 9') (AH) 

and a quantity 

K 6 (g; G) = <U - U(g; G) = [1 + life G)]^ 1 , (A12) 
from which B and / can be simply expressed as 

B Qb (q, q'; G)} = S(q + q') [S ab + IL ab (q; G)} , (A13) 

l ab (q; G) = Y, n ac (g; G) V cb {q; G). (A14) 



1. Schwinger-Dyson equation 

We now derive the Schwinger-Dyson equation for the 
propagator. To this end we have to compute the func- 
tional derivative 5T2[(f>,G]/6G a ^(x,y). Since we only 
consider uniform field configurations, it is simpler to work 
in momentum space, which provides 



5T 2 [<j>,G] A 



5G%(q) 12 A 



S^S ab f ]Tg») 

J a' .. 



A jf Vba (q';G)G%(q + q') A l ab (q;G) ^ b 



1 / A x ^ 

y C.d : p.<7 



xV bd (q';G)G2(q'). 



r c r d V ca (q';G)G%( q + q')x 



(A15) 

Thus, the Schwinger-Dyson equation for the propagator, 
Eqs. (13,14), can be written in terms of vectors and ten- 



's" c,d 
X<t>d[Gdc(q')<t>c] T - 



sors of O(N) as 

[G-\q)] ab [G^ l (q)] ab = ±l N 5 ab ^ Tr (N) G aa (q>) 

+ 377 J q v ba(q';G)G ab ( q + q ') - A \ ab { q -G)4>a4 T b 

3^) I^V ca (q';G)V bd (q';G)G ab (q + q')x 

(A16) 

where (f) a is considered as a column vector and the super- 
script T denotes a transposed quantity. For further con- 
venience, we introduce the function C ab (q) = [G^ 1 (q)] ab 
which is nothing but the 2-point proper vertex. The 
above equations are generalizations to nonzero sources of 
those derived by Mezard and Young 50 within Bray's self- 
consistent screening approximation , which is a zero- 
source version of the 1/N expansion of the 2PI formal- 
ism. 

2. Expansion in free replica sums 

We now expand the different functions entering in 
(A16) in increasing number of free replica sums. The 
general principles of the method are presented in section 
III. The functions G ab (q),C ab (q), U ab (q),V ab (q), I ab (q) 
arc decomposed as in Eq. (22) and their "hat" and "tilde" 
components are expanded as in Eq. (23) and (24). 

For the zeroth-order terms, one easily derives the fol- 
lowing relations: 

G [ \q]<j>a) = C 1 (q;4> a ), 

~[o] -[o] ~[o] -[o] 

G (q;<j) a ,<j) b ) = G {q;<j> a ) C (q\4> a ,<t> b ) G (q;fa), 



(A17) 



V^(q^a) 



[i + n[<%^)]- 



(A18) 



f\ q ;4> a ) = n [ol (g;0 a ) v^ W«), 



f\q;*a,4>b) = -VM(q;<l> a ,<l> b ), 



(A19) 



while from the definition in Eq. (All), the components 
of the polarization can be expressed as 



n[<%^) = A J Tr 0(JV) jc'V^a) 



G l0 \q + q ';<t> a ) + 2G m ( q + q ';<l> ai <t> a ] 



(A20) 
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and 



Ct(q; p) can now be written as 



-jrjy f Tr 0(A r) |g [ 1 (q'l 4>aAb) G [ ° ] {q + q'] <t> a ,4>b) 



(A21) 



A 



G L (q; P r L =q A +m' + 6p+-^ J \G L {q'-p) 



+G L {q': p) + (N- 1) [G T (q'; P) + G T (q'; p) 

V(q'; p){(l- V(q'; p) U(q'; p)) G L (q + q'; p) 



3iV ./,, 



Next, Eqs. (A17)-(A21) are inserted in Eq. (A16) with 
a separate treatment of the diagonal and nondiagonal 
components of the replica matrices. One thus obtains 

two coupled self-consistent equations for C (qi<f> a ) = 

G 1 1 (q;(f> a ) and C [ ] (q;(f> a ,4>b)- 

To proceed further, one must derive the explicit ex- 
pressions for the "bare" or free quantities Cg'(<j;</> a ) = 

G [ ] (q;(j> a ) and Co' {q\<f> a ,<t>b)- From Eq. (A3), one 
finds 



C^iq^a) 



q 2 + m 2 + —\^ a \ 



+G L (q + q';p) 



4A 
3A 



P V(q';p) 2 x 



1 - 2 V(q'; p) U(q'; p)) G L (q'; p) G L (q + q'; p) 



+G L (q'; p)G L (q + q';p) + G L (q' ; p)G L {q + q'; p) 



-4pV(q;p) IL(q;p) 



(A25) 



and 



A 



G T (q;p)- 1 =q 2 + m 2 + 2p+— J {G L (q';p) 



+G L (q';p) + (N-1) [G T {q '; p) + G T (q '; p)) + 



C^ v (q;<i> a> <j> b ) = -A6>»< 



(A22) 



We first consider the 1-replica quantities (i.e. quan- 
tities depending on a single replica field). For explicit 
calculations, one specifies a particular configuration for 

4>a- 



■-^ P^ V(q';p) 2 x 



l-2V(q';p) n(q';p))G L (q';p)G T (q + q';p) 



12pA ! 



(A23) 



+G L (q'; p)G T {q + q; p) + G L (q; p)G T (q + q; p)>, 

(A26) 

where the polarization II(g; p) is given by 



The 1-replica matrices are then diagonal and one intro- 
duces the longitudinal (p = v = \) and transverse com- 
ponents (p = v ^ 1). For instance, 



%q\p) = ^G L {q';p)G L (q + q , ;p)+ 

(N - l)G T (q'; p)G T (q + q'] p) + 2G L (q' ; p)G L (q + q',p) 



G [ ol\q;p) = [G^Uiq^a) =q 2 + m 2 + 6p, 

G^-\q;p) = [Gf-\^ >1 {q;(t> a ) = q 2 +m 2 + 2p. 

(A24) 

For notational simplicity, we drop the superscript [0] 
in the rest of this appendix. The Schwinger-Dyson 
equations for Gh(q\ p)^ 1 = CL{q]p) and Gt(<7;/o) _1 = 



+2(N~l)G T (q';p)G T (q + q\p) 

(A27) 

The expressions of Gi(g;p) and Gr(q;p), which are 2- 
replica quantities evaluated for equal replica fields given 
by Eq. (A23), will be given further below. 

We next consider the 2-replica parts and choose a con- 
figuration of the two replica fields <j> a and </>t> parametrized 
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by 



12pN 



(A28) 



With this choice, the bare 1-replica matrices Co{q',^ a ) 
and Go(q](j> a ) are still diagonal with longitudinal (11) 
and transverse (pp with p > 1) components given by 
Eq. (A24); Co(q;<j>b) and Go(q;<j>b) are not diagonal but 
can be simply diagonalized by applying a transformation 
matrix corresponding to a rotation of angle 9, 



Ro 



For instance, 



cos 9 — sin 9 
sin 9 cos 9 

ljv_ 2 



C (q;(j)b) = Re C (q;p)R-e, 



(A29) 



(A30) 



(diag) ^ (diag) 

with C 011 = C L{q;p') and C 0wi (?;/t/) = 

CoT(q;p') when /i > 1 (and similarly for Go(q;(j>b))- As 
for Co(q; <j> a , (/>&), it is simply given by —A 1, so that 

— ^(diag) ^ (diag) 

Go(?;^,W = -AG (q;p)R g G (q;p)R- 6 . 

(A31) 

The vectorial structure of the correlation functions 
is conserved when fluctuations are taken into account. 
Thus, the previous relations are still true when the bare 
correlation functions are replaced by the "dressed" ones. 
In particular, it is convenient to consider G Re and C Re 
which are related through 

G(q;(j>a,<t>b) Re = 

~(diag) ~ ^(diag) , (A32) 

G [ 9 \q-p)C{q;<S> a ,4>b)Red \q;p'), 

and we introduce the components Gll = (G Re)u, 
Gla = (G Re)i2, Gal = (G Re)2i, Gaa = (G Re)22, 
Gt = {G Re) fj,fj,:fi>2 (and similarly for C), so that 

GLL(q;p,p' \0) = -G L (q;p) G L (q;p')C LL (q;p,p',6) 

G LA (q; P, p', 0) = -G L (q; p) G T (q- p') C LA (q; p, p', 9) 

G AL (q; p, p', 9) = -G T (q; p) G L (q; pf) C AL (q; p, p', 9) 

G AA (q; p, pf, 9) = -G T (q; p) G T (q; pf) C AA (q; p, pf, 9) 

G T (q; p, pf, 8) = -G T (q; p) G T (q; pf) C T (q; p, pf, 8). 

(A33) 

One actually realizes that the description can be fur- 
ther simplified and that only three distinct 2-rcplica 



functions are needed, G^q; p, p' ,£) = G ? ii(g i ; p, p' ,0), 
Gt(q\ p, p' , z) = Gr(q', p, p' ,0), with z = cos 9, as well 
as G A (q; p, p', z) which is such that 



Gla(q;P,p',0) = G AL (q;p',p, -9) = sin6> G A (q; p, p' \z) 
G AA (q; p, p',6) = cos 9 G A (q; p,p',z). 

(A34) 

The Schwingcr-Dyson equations can finally be cast in 
the following compact form: 



GLiq-p)' 1 GLiq-pT 1 5(q + q')-—X 



Q(q + q'\ P, p', z) >G L (q'; p, p', z) - \yfp~p 1 V{q; p)x 



V(q;p')U(q;p,p',z) = Az, 



(A35) 



Gdq-p)- 1 Griq-pT 1 6(q + q')-^X 



Q(q + q'; p, p', z) >G A {q'; p, p', z) = A, 



(A36) 



G T (q; p)- 1 G T (q; pT 1 S(q + q')-^X 



3N 



Q(q + q'; p, p', z) }G T {q'; p, p', z) = A, 



with 



(A37) 



Q(q; p, p',z) = V(q; p) V(q; p'){ U(q; p, p',z)x 



1 — 4/3 Gl (q; p) V(q; p) - V G L (q; pf) V(q; pf) 



+4 Vw 7 G L (q;p, p',z) 



(A38) 
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and 



A 



IL(q;p,p',z) = -^ J ^{N -2 + z 2 )G T {q';p,p',z)x 
G T (q + q'; P, p', z) + G L {q';p, p', z) G L {q + q'; p, p' , z)+ 



(l-z) 



GA(q';p,p',z) G A (q + q';p,p',z)+ 



GaW; p', P, z) G A {q + q'\ p' , p, z) 

(A39) 

The previously introduced functions of one replica field 
G L {q;p), G T (q;p), and R(q;p) (see Eqs.JA25, A26, 
A27)) are simply given by G L {q;p) = G L (q; p, p, 1), 
G T (q;p) = G T (q;p,p, 1), and U(q; p) = Il(q; p, p, 1). 



3. Limit of existence of analytic solutions 

Note that in the limit where the two replicas are equal, 
i.e. p = p' and z = 1 (8 = 0), and provided that 
all functions are regular enough around this limit, one 
obtains from Eqs. (A25- A27) and (A35-A39) a closed 
set of four coupled self-consistent equations for Gi(g; p), 

Gr(<7; p), Gh{q', p) and Gr{q] p). However, as stressed in 
the main text, this breaks down if the field dependence is 
nonanalytic when the two replicas become equal. In the 
present O(N) model, the limit of equal replica fields can 
be approached along several directions: one may for in- 
stance consider the dependence of the 2-replica functions 
on (1 — z) for p and p' fixed and on (y/p— \f(?) for z and 
(y/p + \fp f ) fixed. (Note that the 2-replica functions are 
even in the variable (yfp — y/p') but not in (1 — z).) To 
generalize what has been done in section IV-B for a one- 
component field, we apply d y = d/dy on eqs. (A35-A37) 
with y = z or y = yfp — y/p 7 . So, for instance, one finds 



d v Q(q; p, p' ', z)\ 4 = 4p V{q;p) 



d y G L (q;p,p',z)\,+ 



V{q;pf l-8pG L (q;p)V(q;p) / G T (q + q'- 1 p)x 



under the assumption that the 2-replica functions Gl, 
Ga, Gt, obtained as solutions of the above Schwingcr- 
Dyson equations, are analytic in the limit of equal replica 
fields. It is then easy to realize that Gl = Gt = G, 
Gl = Ga = Gt = G and that, as a result, the operator 
A' ' (<z, q') has a single component corresponding to the 
equation 

/ A [o1 (q, q')d y G(q'; p, p' , z) \ = reg, (A41) 

Jq< 

with "reg" equal to zero if y = yfp — yfp' (due to the 
permutation symmetry between replicas) and "reg" finite 
if y = z. The operator A^(q, q') is found equal to 

AW(q,qi) = G(q;p = 0)- 2 6(q + q')- 



3N 



Q(q + q';0,0,l) 



1 

N 



v(q"; P = oy 



xG(q + q"; P = 0) G(q" + q';p = 0) + O(-). 

(A42) 

In addition, from Eqs. (A38) and (A39), one has 

Q(q; 0, 0, 1) = V(q; p = 0) 2 U(q; p = 0), (A43) 



U(q;p = 0) 



G(q';p = 0) G(q + q';p = 0), 



i' 



(A44) 



with V(q; p = 0) = (1 + U{q: p = 0))" 1 , and 
A 



U(q;p = 0) 



G 



G(q';p = 0)x 



q' 



G(q + q';p = 0) + 2G(q + q';p = 0y j 

(A45) 



It is now easy to check that the operator A^ (q, q'), which 
is associated with the limit of validity of the analytic so- 
lutions to the Schwinger-Dyson equations, coincides with 
the replicon operator derived in Refs. [31,32]. 



d y G T (q';p, p', + <?(-)+ reg, 

(A40) 

where the subscript <fi indicates that the derivatives are 
evaluated for equal replica fields, p = p' and z = 1, and 
"reg" denotes regular terms that involve derivatives of 
1-replica quantities. 

Due to the presence of three different 2-replica func- 
tions, Gl, Ga, Gt, the operator A^(q,q') is now a 
matrix with several distinct components. To avoid too 
much formalism, we derive here the operator in the case 
where p = p' = (and of course z — 1). Generaliza- 
tion is tedious but straightforward. To find the poten- 
tial breakdown of analyticity, we evaluate the operator 



APPENDIX B: STABILITY OPERATOR FOR 
LOCAL SELF-ENERGIES 

In several approximations or models (mean field, Gaus- 
sian variational method, locator and related DMFT ap- 
proximations), the self-energies are taken as purely local 
quantitites. As a result of the Schwinger-Dyson equa- 
tion, Eq. (13), the only momentum dependence in the 
(dressed) propagators comes from their bare counter- 
parts. For instance, in the explicit RSB scheme (con- 
sidering for simplicity a 1-componcnt field), one has 



(Bl) 



15 



By assumption, the self-energies are functions of the lo- 
cal piece of the propagators only, e.g. £[°1(0, 4>') only 
depends on G^(x = 0; = J q G^(q; <t>, <j/). 

From Eq. (37), the stability operator A^°\q, q'\ <p, cf) = 
Ai?(<7, q') is then given by 



A R (q,q')=G(q;<t>)- 2 6(q+q')- 



SG(x = O;0, V) 



(B3) 



where the last term has no momentum dependence and 
is evaluated for <j>' = (p. 

Finding the eigenvalues of this stability operator is now 
an easy task. There is in fact a single eigenvalue Xr which 
is solution of 



1 = 2 

SG 



G{q-Af 



(B4) 



tJ q l-\ R G{q;<j>) 2 
with the associated eigenfunction ^S R (q;<p) oc 



G(q;4>) 2 (l - X R G(q; 0) 2 )" 1 (above and in what fol- 
lows, we drop the superscript [0]). The operator 
becomes marginal when X R = 0, which characterizes 
a generalized Almeida-Thouless line 1 ' and happens 



(B2) whenever 1 = SE/SG^ J G{q; 



The inverse operator XR(q>q') is obtained from 
Eq. (B3) as 



XR(q,q') = G(q)' 



S(q + q') + 



SG \<t> 



G{q'f 



(B5) 

and the spin-glass susceptibility, \SG = J q f q , Xfl(g, q'), 
is simply given by 



XSG 



I G{q? 



(B6) 



In this case, which encompasses most approximations 
through which spontaneous RSB has been found, the 
spin-glass susceptibility diverges as soon as the replicon 
eigenvalue goes to zero. 
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